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Abstract 


We investigate the stabilization problem of a cascade of a fractional ordi- 
nary differential equation (RODE) and a fractional diffusion (FD) equation, 
where the interconnections are of Neumann type. We exploit the PDE back- 
stepping method as a powerful tool for designing a controller to show the 
Mittag—Leffler stability of the FD-FODE cascade. Finally, numerical simu- 
lations are presented to verify the results. 
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1 Introduction 


One of the most useful approach to obtain a boundary controller is the par- 
tial differential equation (PDE) backstepping method, which was initially 
utilized in [5, 6] for spatial discretization case and has since been expanded 
for continuous case in [17] for many applications, like fluid flows [1, 2]. The 
PDE backstepping approach, which is an integral operator with a known 


*Corresponding author 

Received 1 January 2020; revised 28 February 2020; accepted 5 March 2020 

Shadi Amiri 

Department of Applied Mathematics, Faculty of Mathematical Sciences, University of 
Guilan, Rasht, Iran. e-mail: amiri.shadi63@gmail.com 


Mohammad Keyanpour 

Faculty of Mathematical Sciences, and Center of Excellence for Mathematical Modelling, 
Optimization and Combinational Computing (MMOCC), University of Guilan, Rasht, 
Tran. e-mail: kianpour@guilan.ac.ir 


+ This article was suggested by the Scientific Committee of the Third Iranian Sem- 


inar on Control and Optimization for publication in IJNAO, which was accepted after 
independent review. 


177 


178 Sh. Amiri and M. Keyanpour 


and continuous kernel function, has been used to analyze some problems of 
boundary stabilization of integer order unstable heat system with boundary 
control law. For example, this method was used in [21] for boundary feed- 
back stabilization of heat equation. Also, this approach was applied in [26] 
to investigate the boundary stabilization of a class of linear parabolic partial 
integro-differential equation. 

Since real systems in our world are complex and display memory and 
genetic characteristics, they can be well characterized by fractional order’s 
notions; see [29, 22]. On the other hand, it is confirmed that fractional cal- 
culus [29] is very effective in modeling and analysis of distributed parameter 
systems [8]. In the last decades, researchers have interested in investigat- 
ing the stability of the fractional-order systems. In [20], the stability of the 
fractional-order linear system subject to input saturation has been discussed. 
Since most of the systems in the real world are nonlinear, the stability prob- 
lem of the nonlinear systems has become an attractive issue to study. For 
example, in [7, 10, 18], the Mittag—Leffler stability of the fractional-order 
nonlinear systems for 0 < a < 1 was addressed. Authors in [13] focused 
on stability analysis of a class of fractional-order nonlinear systems with 
0O<a<2. 

Applying the backstepping method in designing a controller for fractional 
ordinary differential equations (FODEs), was first used in [11]. For instance, 
one can refer to the work [9] in which the adaptive fractional-order backstep- 
ping was used to design an adaptive feedback control law that Mittag—Leffler 
stabilize the commensurate fractional-order nonlinear systems. On the other 
hand, some researches have been dedicated to the fractional-order PDE sys- 
tems, where the backstepping method is applied to design a controller to 
solve the stabilization problem of the mentioned systems. However, com- 
paring to the ordinary fractional-order systems, the stabilization problem 
of fractional-order PDE systems has been less investigated. For example, in 
[23], the stability problem of one-dimensional wave equation was discussed via 
boundary fractional derivative control, and in [12], the backstepping method 
was applied to investigate the stabilization problem of the fractional diffu- 
sion (FD) system, governed by the FD equation consisting the diffusion term, 
with Dirichlet or Neumann condition. In [8], the boundary feedback control 
problem of the FD system with mixed or Robin boundary control was ad- 
dressed via the backstepping method. In [30], the backstepping method was 
used for the stability problem of a class of unstable time fractional diffusion 
equation with the Dirichlet and Neumann boundary controls. 

According to the previous paragraphs, systems described by ordinary dif- 
ferential equations (ODEs) and also systems modeled by PDEs are common 
in control engineering, and many works have been dedicated to the theory 
of them. Recently, the coupled systems have become one of the interesting 
areas of study. Examples of these systems are provided in control problems 
of electromagnetic coupling, mechanical coupling, and chemical reaction cou- 
pling [28]. In the last decades, the stabilization problem for coupled systems 
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has become challenging areas. The cascade structure for the heat PDE with 
an ODE and also the cascade structure for a wave PDE with an ODE, when 
the interconnection is of Dirichlet type, were discussed in [15] and [16], re- 
spectively. Next, these results for the stability analysis of the PDE-ODE 
cascade when the interconnection is of Neumann type have been extended by 
Susto and Krstic [27]. In [28], the stabilization problem for a new cascade of 
PDE-ODE was studied. 

To the best of our knowledge, the stabilization problem and also con- 
troller design for a cascade of an FD equation and an FODE have been less 
investigated. In this paper, we consider a cascade of an FODE and an FD 
equation, and we use an invertible integral transformation to transfer the 
original system to a Mittag—Leffler stable target system. Finally, we present 
a numerical example for verifying the results. 

Notation: L?(0,D) represents the usual Lebesgue integrable functions. 
Let u(.,t) € L?(0, D); then we define 


D 
||u(., t)|| = (/ Hepa) 


Also, we denote a symmetric negative definite matrix A € R"*” by A ~ 0. 


1 


2 


2 Preliminaries 


Definition 1. [4] The Caputo fractional-order derivative is defined by 


ener ee * Xx \(r) 
PO =e [ (f—r)ati-n (n-—l<a<n) (1) 


where a is the order of fractional derivative and the gamma function I is 
defined as (7) = Jg° t7~*e~*dt. 
Definition 2. [14] The Mittag—Leffler function is defined as 


ok 


Bal?) = sa ans 
=) T(ka +1) 
where 0 < a < 1. The Mittag—Leffler function with two parameters is given 
by 


ok 


Ea,p(2) = > Tika A) 


Definition 3. [19, 18](Mittag—Leffler stability) The solution of 


180 Sh. Amiri and M. Keyanpour 
fy Deu(t) = f(t,u), 
is said to be Mittag—Leffler stable if 
ayy 
Ilu(t) |] < (mluto)]Ba(-AlE — to)*)) 
where ¢o is the initial value of time, a € (0,1), A > 0, b > 0, m(0) = 0, and 


m(u) is nonnegative and meets locally Lipschitz condition on u € B Cc R” 
with the Lipschitz constant mo. 


Lemma 1. [3] Suppose that u : [0,co) > R is a continuous and differ- 
entiable function. For any time t > to > 0, one can readily show that 


1 
~¢ Deu? (t) < u(t)g Diu(t), O<a <1. 


Lemma 2. [13] Assume that X : [0,co) > R” is a vector of differentiable 
function. If a continuous function V : [tp,00) x R” > R satisfies 


EDFV (ts XO) < —yV (t, X(t), 
then 
V(t, X(B) < V(to, X (to)) Baa(—7(t — to)*), 


where 0 < a < 1 and 7 is a positive constant. 


3 Problem statement and analysis 


Consider the cascade of fractional diffusion (FD) and a fractional-order ordi- 
nary differential equation (FODE) with Caputo derivative as follows: 


oD? X(t) = AX(t) + Buz(0,t), (2) 
(DPC) = tee(a2), (3) 
u(0,t) = 0, (4) 
u(D,t) = U(), (5) 


where X(t) € R” and u(z, t) are the state of the FODE and FD, respectively, 
and U(t) is the control input. Note that t > 0 and a € [0, D] in which D > 0 
is the length of the PDE domain. The aim is to Mittag—Leffler stabilize the 
system (2)—(5). 

The PDE backstepping, introduced by Krstic, is the most effective ap- 
proach for boundary controller designing for the PDE systems. In this 
method, we use an invertible integral transformation (X,u) — (X,w) to 
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convert the cascade of an FD and an FODE (2)-—(5) into the following target 
system: 


6 D¢X(t) = (A+ BK)X(t) + Bu,(0,t), (6) 
0 Drw(z,t) = Wee (2, t), (7) 
w(0,t) = 0, (8) 
w(D,t) = 0. (9) 


The control gain K is chosen such that the Mittag—Leffler stability of the tar- 
get system is guaranteed. Since the transformation is invertible, the Mittag— 
Leffler stabilization of the original closed-loop system will be derived. 


4 Designing the state feedback controller 


We consider the following transformation: 


w(z,t) = ula, t) — i q(x, y)uly, t)dy — y(2) X(t), (10) 


where the unknown functions q(x, y) and (a) should be determined to con- 
vert system (2)—(5) into the target system (6)—(9). First, we derive the ker- 
nels, and next, we prove that the target system (6)—(9) is Mittag—Leffler 
stable. 


To determine the unknown functions, we need the first two derivatives of 
w(x, t) with respect to x that are given by 


te(t,t) = u(r, t) — q(a,2)u(e, t) — | ‘gaauy tay FORO, GD 


Wea (#,t) = Une (x,t) — (g(a, a))'ula, t) — q(a, a) ta (aw, t) — dew, aula, t) 


= [ dua(x, y)u(y, t)dy — v(x) X(t), 
(12) 


and we take the Caputo fractional derivative of w(z,t) respect to ¢: 


182 Sh. Amiri and M. Keyanpour 
Hi 
© Dew( x,t) = ©Deu(x,t) — | ole, 9) Deu(y, t)dy — (x) DEX (t) 
0 


= tna) = f° a(eudiny(vst)dy 1 )(AX(W + Bus(0,0) 
= tne (t,t) — a(t, 2) tte, ) + fale, 0) — (2) Blese(0,1) 


+ Gy (x, x)u(a, t) — [ duy (2; y)uly, t)dy — y(@)AX(#). (13) 


Now we evaluate the backstepping transformation (10), and (11) in x = 0 
and we also subtract the right-hand side of (12) from the right-hand side of 
(13), then 


w(0,t) = —7(0)X (t), 
We (0, t) = uz(0,t) — 7/(0) X(t), 
0 Drw(x,t) — wee(x, t) = 2(q/(x, x))u(z, t) + (g(x, 0) — 7(x) Bluz(0, t) 
+ [ [dex (2,4) — dyy (x,y) u(y, t)dy 
+ [y"(2) — y(a) A] X(t), 


in which we have used the Dirichlet boundary condition u(0,t) = 0. A 
sufficient condition for equations (7)—(9) to be held is that the unknown 
functions (a) and q(z, y) satisfying an ODE of second order and a hyperbolic 
PDE of second order that come, respectively, 


y"(x) = Ay(2) (14 
7(0) = 0, (15 
7'(0) = K, (16 
and 
ua (X,¥) = dyy(2,¥), (17 
q(x, x) =0, (18 
q(z,0) = y(2)B (19 


According to [27], the solution of (14)—(16) is 


V(#) = KM(a) (20) 


and the solution to (17)—(19) is 
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q(x, y) = s(z — y) 


in which we have used the functions s(.) and M(.) for simplifying of notation 
in the proof. 


By [27], the backstepping transformation is invertible, and the inverse 
change of variables is as follows: 


u(x,t) = w(a,t) + - r(x, y)w(y, t)dy + A(x) X(t), (22) 


where the kernel functions r(z,y) € R and A(z) € R” are determined as 
follows: 


A"(z) = (A+ BK)X(z), (23) 
(0) = 0, (24 
(0) = K, (25 
and 
Tax (x, y) = lyy (a, Y), (26 
r(x,x) = 0, (27) 
r(x, 0) = (x) B (28 
Therefore, the solutions of (23)-(25) and (26)—(28) are, respectively, 
Ma) = KN(2), 
and 
r(z,y) = n(a,y), 
where 
i. 
w@=lorell 9 Fil, 
n(f) = KN(E)B. 
Also, based on relations (10) and (22), we can write: 
we(2) = us(e) — [sale —wudy-KM'@XO), 9) 
Ux (Z) = We(x) + a Ny(x — y)w(y)dy + KN'(x) X(t). (30) 


According to [27], from (29) and (30), we have 
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||we||? < ar ||uel|? + ag ||ull? + a3||X|]?, (31) 
[[ue ||? < Bal|eve||* + Bal|evl|? + B3llX|/?. (32) 


Now we prove that the target system (2)—(5) is Mittag—Leffler stable. Then 
with the help of the backstepping method and using an invertible transfor- 
mation, we obtain the Mittag—Leffler stability of the original system (6)— 
(9). Before stating the theorem, we need to consider the following lemma. 
Lemma 3. There exist positive constants j4; and p12 such that 


|[etae|I? + IX]? < Ha (llwe ll? + |X 117), (33) 
and 
||w2 (a, 0)||? + |X (0) |? < p2((|e(x, 0)||? + |X (0) |I?). (34) 


Proof. By using the Poincaré inequality for ||w||?, we can write (32) in the 
following form: 


Jus ||? < max {61 + 4D* Bo, 83} (\lwell? + |X|”). (35) 


Let 
p = max { 8; + 4D” Bo, 83}, 


then, it is clear that: 
[|x |? + |X]? < pa(|fwol|? + X17), 


in which yw; = p; + 1. With the help of (31) and in a similar manner, (34) is 
obtained. 


Also, we consider the following assumption throughout of the paper: 


(H1) We assume that system (2) is controllable. 


Theorem 1. Consider the closed-loop system (2)—(5) with the control law: 


n Jo 


: feats o [in tlomn, 
U(t) = K [On In] {¢ E Jxo+f e Fa Bu(ysthay} (36) 


Assume that there exist positive constants d and / and also a symmetric pos- 
itive definite matrix P, such that the control gain K satisfies in the following 
inequality: 


P(A+BK)+(A+BK)?PPB 0 0 
BTp -d 0 0 
Q= 0 0 —B 0 <0. (37) 


0 0 0 -6 
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Also, suppose that u,(2,0) is square integrable in x for any initial condition. 
If w,(.,t) is continuously differentiable on t € [0,00], then the closed-loop 
system under the control law (36) is Mittag—Leffler stable under the norm 


(XI? + fo u2 (a, t)dx)2. 


Proof. We consider the following Lyapunov function: 
7 c 2,4 2 
V(t) = XTPX + Spl? + $level? (38) 


in which c > 0 and d > 0 will be chosen later. Also, ||w(., ¢)||? and ||w.(.,¢)||? 
represent the notation for Lz norms fe w(x, t)dx and he w2 (x, t)dax, respec- 
tively. By taking the Caputo derivative of V with respect to t, we have 
c d 
DEV =§DE(XTPX + Shell? + 5 Ihe”) 


a d Qa 
o Dewi? +5 6 De llwell?. (39) 


Cc 
Sipe x Px) 5 0 


2 


According to Lemma 4 in [8], w(z,t) is continuous and differentiable on 
t € [0,co). On the other hand, by the assumption, w,(.,¢) is continuously 
differentiable, so by Lemma 1, we have 


D D D 
gpe [ w* (2, t)dx = / 6 Dew? (a, t)dx < 2 | w(a, t)f Dew(a, t)dx 
0 0 0 


=2 i; w(x, t)Wan(x, t)dx 
0 
= —2\well?, 
in which we have used integration by parts, so: 
o De llwll? < —2||well- (40) 


To compute the Caputo derivative of ||w,.||?, we multiply wz(2,¢t) by (7) and 
integrating from 0 to D. Then 


D D 
[ waS Dew(e,de= [ w2,(e,t)de, (41) 
0 0 


and applying integration by parts to the left side of (41), we have 


D D 
| Wee 6 Dow(z, t)de = 0 — i We 5 Dew,(2, t)dz. (42) 
0 0 

Because of the Dirichlet boundary condition w(0,t) = 0 and w(D,t) = 0 and 


based on the Caputo time fractional derivative’s definition in [25], we have 
o Dew(0,t) = 7 Dew(D,t) = F De0 = 0, (for all t € [0,00)). By considering 
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(41) and (42), we have 


»D D 
/ We(z, t)S De w(a, t)dx = -{ we, (x, t)dx. 
0 0 
Now we can evaluate the Caputo derivative of ||wz||? respect to t as follows: 
D D 
C De ||wel|? = gps | w2(a,t)dr = CD? w2 (a, t)dr 
0 0 
D 
< 2 | we(a,t) De we(2x,t)dx = —2||weell?. (43) 
0 
On the other hand, applying Lemma 1 and using (6), it is easy to see that 
CpnayT T T 
C p2xXT PX <XT(t) (PU + Bk) +(A+ BK) P)X(t) 
+ BT PXw;,(0,t) + X? PBw,(0, t). (44) 
Here, the goal is to establish § D?V(t) < —yV/(t), for some positive constant 
+ to prove the Mittag—Leffler stability of the target system (6)—(9). Therefore, 
we replace relations (40)—(44) into (39), and we have 
CO DeV <X7(P(A+ BK) +(A+BK)'P)X 
+ BY’PXw,(0,t) + X7 PBw,(0,t) — ¢llwa||? — d||weell?. 


It can be shown by Agmon’s inequality that for the system (6)—(9) the fol- 


lowing inequality holds: 
1+ D 
= ||teal|? < —5—|lwe|? — w; (0). (45) 


x 
Hence 


ODeV < X7(P(A+ BK) +(A+BK)'P)X 
1+D 


+ B’PXw,(0,t) + X? PBw,(0,t) — (c— d——)|lwell” — dw? (0, t). 
(46) 
Now, we assume that 
1+D 
f= So), 47 
c D > (47) 


By using Poincaré inequality, we can rewrite (46) in the following form: 
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OD2V <X7T(P(A+ BK)+(A+ BK)'P)X 
1 1+D 


BTPX f)4+X7 PB t wl” 
+ we(0,t) + we (0,4) — apa le pel 
1 1+D es 
-—d We||,? —dw2(0, t). 48 
Tp ap? DH) llwell” —adwe (0, t) (48) 
We can rewrite inequality (48) as follows: 
X |° [P(A+BK)+(A+BK)™P PB O 0 xX 
Cna w, (0, t) BE P -d 0 0 w,, (0, t) 
DEV < 
oer 7 [wll 0 0 =2°0 eal 
||we| 0 0 0 —B} | [well 
49 


Assumption (37) and relation (49) imply that fD@V <0. Since 2 < 0 and 
by (47), we should have 


Ge (c iF") >0, (50) 


~ 14 4D? D 


and we can conclude that 


FeV =o, 


Let Ao := Amin(—2). Then »» > 0, and for any nonzero values of X(é), 
w(a,t), and w,(z,t), we have 


0 DEV < —o( (|X|? + [lve (0, #)||? + [lew]? + [lee ll?) 
< —o (|X|? + |lewll? + [we l?). (51) 


Since P is a positive definite matrix, then 
T i 2,4 2 ove 2,4 2 
V = XT™PX + Heol? + Seve < Amae(P)IXI? + Stell? + §lrvell 
< max { Anns(P)s $5 (CTP + ll? + lta) 
so we can rewrite inequality (51) in the following form: 


Dev <-yW 


Xo 
max{\max(P),$, 2} 


in which y = . It follows from Lemma 2 that 
V <V(0)Fei(—yt*). (52) 


Therefore, we can write 
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2Amin(P)|X ||? + dllwel|? < 2X7 PX + ellw||? + dllwel|? < 2V(0)Ea(—7#t*) 


and 


2 
min {2Amin(P), d} ‘ 


IX? + llwell? < (0) Bo(—9#t*). (53) 


On the other hand, using the Poincaré inequality in the Lyapunov function 
(38), we have 


c d d 
V= XPPX + sllwl? + sllwell? < max { Snax(P) 200? + $b (XI? + wel) (54) 


Therefore, by (54), we can write (53) in the following form: 


2k 


x|l? all? < 
I I + Iw ] = min {2Amin(P), d} 


(Xl? + llwe(#, 0)|/?)Ba(—7t*), (55) 


in which 


R= max {Amax(P), 2cD? + 5} ’ 


By taking m(w) in Definition 3 as 

m((||X (0) |}, [lve (2, 0)II)) = saga ay (IX (0) |? + [lw2(2, 0)|I?), 

it is clear that m/((||X(0)||, ||we(z,0)||)) is locally Lipschitz in ||X(0)|| and 
||w.. (ax, 0)||?. On the other hand, m/((||X(0)||, |/w(x, 0)|])) > 0 for ||X(0)|| 4 0 
or ||wz(x,0)|| 4 0 and it is zero if and only if || (0)|| and ||wa(z, 0)|| are zero. 
Then by (55) and Definition 3, the target system (6)—(9) is Mittag—Leffler 
stable under the norm ||X(t)||+||w.(a, t)||?. 

Now, by Lemma 3 and relation (33), we can write (55) in the following form: 


2A 


Fale all? < 
[AI + Metall Ses {Dm (P) 


(|X (0)? + |]w2(x, 0) ||?) Za(—7t"), 


then from (34), we have 


2Ru1 be 
min {2\min(P), d} 


|X|? + lluell? < (|X(0)||? + ||ue (x, 0) ||?) Za(—7t*), 


(56) 


which guarantees that the original system (2)—(5) is Mittag—Lefller stable. 


5 Numerical simulation 


In this section, we present a numerical example to verify our theoretical re- 
sults. In this example, we discuss the stability Mittag—Leffler results related 
to system (2)—(5). 
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Example 1. We consider the following unstable system: 


0 Dr X(t) = X(t) + uz(0,¢), (57) 
DP il, ) = tpg Gt) (58) 
u(0,t) = 0, (59) 
u(1,t) = U(t) (60) 


Comparing with the coupled system (2)-—(5), it is clear A = 1, B = 1, and 
D = 1. We assume that a = 0.75 and u(#,0) = 0 and X(0) = 1 as initial 
conditions. Also, U(t) is determined by relation (36) in the following form: 


U(t)=K sin) x( Pp] sinh(1 — y)u(y)dy. (61) 


In order to show the stability of the closed-loop system (57)—(60) under the 
control law (61), we select c = 7 and d = 1 in the Lyapunov function (38) to 
satisfy (47) and we also choose P = 30. Then, we obtain the feedback gain 
K, with the help of CVX 1.2.1, to satisfy the stability condition: 


P(A+BK)+(A+BK)’P PB 0 0 


BP —-d 0 0 er, 
0) 0 -6 0 : 
0 0 O -£ 
We get 
kK = —16.0150. (62) 


To figure out the state variable X(t), u(x,t) of the cascaded FODE-FD 
system (2)—-(5), we have used the finite-difference approximation method de- 
scribed in [24] to descretize the spatial solution domain [0, D] into finite 
numbers of Q + 1 subintervals and the time interval [0, T] into TM +1 grid 
points for some positive integers Q and M, that is, this finite difference al- 

1 


gorithm estimates the system (2)-(5) with the special stepsize h = G and 


k= 4G for x and t, respectively. 


Now, we set the discretization parameters T = 20, Q = 50, and M = 64. 

With these discretization parameters and considering (62), we can see 
from Figure 1 that the state of the coupled FODE-FD system, that is, X(t) 
and u(x,t), converges to zero for all x € [0, D] and implies that the closed 
system is Mittag—Leffler stable. 

On the other hand, if one selects K = 2 and parameters c, d as before, 
then with the help of CVX 1.2.1, we found out that there is no P that sat- 
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isfies inequality (37) , so the stability condition (37) is not hold for K = 2. 
With the same discretization parameters, we obtain the trajectory of state 
X(t) and u(z,t) of the coupled system (57)—(60). Figure 2 shows that the 
cascade (57)—(60) with the boundary control law (61) is unstable. 


Therefore, one can result that the Mittag—Leffler stablity of the state 
of the closed-loop coupled system (57)—(60) is guaranteed by choosing K 
satisfying the stability condition (37). 


u(x,t) 


02+ 7 : : _— — 20 


0 5 10 15 20 0 0 Time 
Time 


Figure 1: Evolution of state in Example 1. for K = —16.0150. 


. x10 


0.5 


0 5 10 15 20 
Time 


Figure 2: Evolution of state in Example 1 for kK = 2. 
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6 Conclusions 


In this article, we applied the backstepping method in which an invertible in- 
tegral transformation is used and we designed a Dirichlet boundary feedback 
control to guarantee the Mittag—Leffler stability of the cascaded FD-FODE 
system. We presented a numerical example to confirm the obtained results. 
Since the stability analysis of the FODE-FD coupled system is less addressed 
in the literature, this paper is a beginning for the development of the sta- 
bilization problem of the cascade of a fractional-order ordinary differential 
equation and a fractional diffusion equation. 


References 


1. Aamo, O.M., Smyshlyaev, A., and Krsti¢é, M. Boundary control of the 
linearized Ginzburg-Landau model of vortex shedding, SIAM J. Control 
Optim. 43 (2005), 1953-1971. 


2. Aamo, O.M., Smyshlyaev, A., Krstic, M. and Foss, B.A. Stabilization of 
a Ginzburg-Landau model of vortex shedding by output feedback bound- 
ary control, Decision and Control, 2004. CDC. 43rd IEEE Conference, 3 
(2004), 2409-2416. 


3. Aguila-Camacho, N., Duarte-Mermoud, M. A. and Gallegos, J.A. Lya- 
punov functions for fractional order systems, Commun. Nonlinear Sci. 
Numer. Simul. 19 (2014), 2951-2957. 


4. Ahn, H.-S., Chen, Y.Q. and Podlubny, I. Robust stability test of a class 
of linear time-invariant interval fractional-order system using Lyapunov 
inequality, Appl. Math. Comput. 187 (2007), 27-34. 


5. Balogh, A. and Krstic, M. Infinite dimensional backstepping-style feedback 
transformations for a heat equation with an arbitrary level of instability, 
Eur. J. Control. 8(2) 2002, 165-175. 


6. Balogh, A. and Krstic, M. Stability of partial difference equations govern- 
ing control gains in infinite-dimensional backstepping, Systems Control 
Lett. 51 (2004), 151-164. 


7. Chen, J., Zeng, Z. and Jiang, P. Global Mittag-Leffler stability and syn- 
chronization of memristor-based fractional-order neural networks Neural 
Networks, 51, (2014), 1-8. 


8. Chen, J., Zhuang, B., Chen, Y.Q. and Cui, B. Backstepping-based bound- 
ary feedback control for a fractional reaction diffusion system with mixed 
or Robin boundary conditions, IET Control Theory Appl. 11(17) (2017), 
2964-2976. 


192 Sh. Amiri and M. Keyanpour 


9. Ding, D., Qi, D., Meng, Y. and Xu, L. Adaptive Mittag-Leffler stabilization 
of commensurate fractional-order nonlinear systems Decision and Control 


(CDC), IEEE 53rd Annual Conference, (2014), 6920-6926. 


10. Ding, D., Qi, D. and Wang, Q. Non-linear Mittag-Leffler stabilisation of 
commensurate fractional-order non-linear systems, IET Control Theory 
Appl. 9 (2015), 681-690. 


11. Efe, M.O. Application of backstepping control technique to fractional or- 
der dynamic systems, Fractional dynamics and control, 33-47, Springer, 
New York, 2012. 


12. Ge, F., Chen, Y.Q. and Kou, C. Boundary feedback stabilisation for the 
time fractional-order anomalous diffusion system, IET Control Theory 


Appl. 10 (2016), 1250-1257. 


13. Huang, S. and Wang, B. Stability and stabilization of a class of fractional- 
order nonlinear systems for 0 < a < 2, Nonlinear Dynam. 88 (2017), 
973-984. 


14. Huang, S., Zhang, R. and Chen, D. Stability of nonlinear fractional-order 
time varying systems, J. Comput.Nonlin. Dyn. 11(3) (2016), 9 pp. 


15. Krstic, M. Compensating actuator and sensor dynamics governed by dif- 
fusion PDEs, Systems Control Lett. 58 (2009), 372-377. 


16. Krstic, M. Compensating a string PDE in the actuation or sensing path 
of an unstable ODE, IEEE Trans. Automat. Control 54 (2009), 1362-1368. 


17. Krstic, M. and Smyshlyaev, A. Boundary control of PDEs: A course on 
backstepping designs SIAM, 2008. 


18. Li, Y., Chen, Y.Q. and Podlubny, I. Mittag-Leffler stability of fractional 
order nonlinear dynamic systems, Automatica J. IFAC 45 (2009), 1965- 
1969. 


19. Li, Y., and Chen, Y.Q. and Podlubny, I. Stability of fractional-order non- 
linear dynamic systems: Lyapunov direct method and generalized Mittag- 
Leffler stability, Comput. Math. Appl. 59 (2010), 1810-1821. 


20. Lim, Y.-H., Oh, K.-K. and Ahn, H.-S. Stability and stabilization of 
fractional-order linear systems subject to input saturation, IEEE Trans. 
Automat. Control 58 (2013), 1062-1067. 


21. Liu, W. Boundary feedback stabilization of an unstable heat equation, 
SIAM J. Control Optim. 42 (2003), 1033-1043. 


22. Mandelbrot, B.B. The fractal geometry of nature, Schriftenreihe ftir den 
Referenten. [Series for the Referee] W. H. Freeman and Co., San Francisco, 
Calif., 1982. 


On the stabilization of a coupled fractional ordinary and partial ... 193 


23. Mbodje, B. and Montseny, G. Boundary fractional derivative control of 
the wave equation, IEEE Trans. Automat. Control 40 (1995), 378-382. 


24. Murio, D.A. Implicit finite difference approximation for time fractional 
diffusion equations, Comput. Math. Appl. 56 (2008), 1138-1145. 


25. Podlubny, I. Fractional differential equations: An introduction to frac- 
tional derivatives, fractional differential equations, to methods of their so- 
lution and some of their applications, Elsevier, 1998. 


26. Smyshlyaev, A. and Krstic, M. Closed-form boundary state feedbacks for 
a class of 1-D partial integro-differential equations, IEEE Trans. Automat. 
Control 49 (2004), 2185-2202. 


27. Susto, G.A. and Krstic, M. Control of PDE-ODE cascades with Neumann 
interconnections, J. Franklin I. 347(1) 2010, 284-314. 


28. Tang, S. and Xie, C. Stabilization for a coupled PDE-ODE control system, 
J. Franklin Inst. 348 (2011), 2142-2155. 


29. Torvik, P.J. and Bagley, R.L. On the appearance of the fractional deriva- 
tive in the behavior of real materials, J. Appl. Mech. 51(2) 1984, 294-298. 


30. Zhou, H.-C. and Guo, B.-Z. Boundary feedback stabilization for an un- 
stable time fractional reaction diffusion equation, SIAM J. Control Optim. 
56 (2018), 75-101. 


